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A small subset of normal modes mimics the properties of dynamical 
heterogeneity in a model supercooled liquid 

Glen M. Hocky 1 and David R. Reichman 1 ^ 

Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, 
USA 

In this work, we study the nature of transitions between inherent structures of a two-dimensional model 
supercooled liquid. We demonstrate that these transitions occur predominately along a small number of 
directions on the energy landscape. Moreover, we show that the number of such directions decreases as the 
temperature of the liquid is decreased in the mildly supercooled regime, in concert with earlier studies on an 
athcrmal jamming system. We show that this decrease happens in parallel with a change in character of the 
transitions as dynamics in the system becomes more heterogeneous and localized. We investigate the origin 
of these trends, which suggests interesting connections between jamming and thermal glassy phenomena. 



I. INTRODUCTION 

Understanding the formation of glass from a slowly 
cooled liquid remains one of the great unsolved prob- 
lems in condensed matter science i - — The lack of an ob- 
vious change in symmetry and the fact that the transfor- 
mation does not occur at a well-defined thermodynamic 
critical point renders the problem of vitrification more 
difficult to describe than crystallization. Despite this 
challenge, many aspects of glass formation appear to be 
generic and not tied to the specific details of the liquid 
under investigation. In particular, the phenomenology 
of dynamical heterogeneity j 4 - which includes growing dy- 
namical length scales, violations of the Stokes-Einstein 
relation, and prominent non-Gaussian displacements in 
the tails of real-space distribution functions, provides a 
framework upon which theories of the glass transition 
must be based£ 

A problem that bares some similarity with the stan- 
dard laboratory glass transition is the jamming transition 
of hard spheres^ - — Here, slow compression of the sys- 
tem may lead to the formation of a disordered solid. This 
problem is simpler than the vitrification of typical liq- 
uids in the sense that a single control variable, the pack- 
ing fraction, unambiguously tunes the transition upon 
approach to an amorphous configuration with a maxi- 
mally allowed density. While the relationship between 
the standard glass transition and the jamming transi- 
tion is currently vigorously debated^ - — it is clear that 
many aspects of the behavior of supercooled liquids and 
suspensions close to the jamming transition share impor- 
tant similarities. In particular, jamming systems display 
the hallmarks of dynamical heterogeneity first exposed in 
finite temperature fluids^ - — 

One approach that has been useful in both the study 
of supercooled liquids and the jamming transition in- 
vokes the notion of soft modest - — Soft modes repre- 
sent a subset of the low frequency harmonic or quasi- 
harmonic displacements of particles around metastable 
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configurations ! 18 i 19 Within traditional mean-field theo- 
ries of glassy liquids, soft modes characterize the mo- 
tion that involves relaxation of groups of particles when 
marginally metastable states first appeari 3 ^ - — Jam- 
ming systems provide a concrete case where such modes 
play a prominent role. In particular, numerical simula- 
tions provide clear indications of diverging length scale 
connected to the jamming transition associated with soft 
modest In addition, soft modes appear to play a criti- 
cal role in defining the mechanical stability of athermal 
packings close to jamming ! 23 ' 24 Currently the evidence 
for a primary role played by soft modes in the jamming 
transition is more extensive than it is for vitrification in 
thermal systems. 

Despite the aforementioned lack of clarity of the roll 
played by soft modes in supercooled liquids, interest- 
ing correlations between dynamics and soft modes in 
such systems have been uncovered via computer sim- 
ulations. In a pioneering set of studies, Harrowell 
and coworkers defined the "isoconfigurational ensem- 
ble" as a means of quantitatively assessing the influence 
of structure on subsequent dynamics £5r— These stud- 
ies revealed that regions that were quantifiably "softer" 
than average were more likely to be involved in large 
amplitude dynamically heterogeneous motion when av- 
eraged over many independent simulations^ - — Later 
it was demonstrated that the real space properties of 
the inherent structure normal modes correlate closely 
with dynamical heterogeneities and irreversible configu- 
rational rearrangements^ - — Such connections have also 
been demonstrated in experiments performed on colloidal 
suspensions^ - — It remains unclear in the examples dis- 
cussed above if the harmonic soft modes are an active 
player in the dynamics or if they are spectators whose 
positions correlate with regions of dynamical activity but 
which do not influence particle relaxation. 

With an appropriate definition of normal modes 
emerging from the free energy landscape of hard spheres, 
Brito and Wyart have shown that dynamical heterogene- 
ity in jamming systems also correlates with low frequency 
modes.— They have demonstrated that relaxation events 
in simulated hard spheres proceed along a small number 
of mode directions. Further, they have explicitly demon- 
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FIG. 1. Radial distribution functions for the 2D 65:35 Kob- 
Andersen Lennard- Jones system at p — 1.2. Temperatures 
are T = {5.0, 3.0, 2.0, 1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.45} from bot- 
tom to top. Each g(r) is shifted up by 0.5 from that of the 
preceding temperature for clarity. 

strated that the number of such directions systematically 
decreases as the jamming threshold is approachedi 24 ' 36 
The purpose of this work is to explore the possibility 
that a similar effect may occur in thermal systems as 
temperature is lowered, even in the absence of an analog 
of the jamming threshold where the correlation length 
associated with soft modes diverges. 

The paper is organized as follows. In Section|TT]we will 
present the details of the two-dimensional glass forming 
liquid that we study here. In Section Mil we will review 
the notion of inherent structures (IS) and IS trajecto- 
ries, and present relevant properties of the IS trajecto- 
ries for our model. In Section HVl we will discuss normal 
modes (NM), the NM properties of our system, and the 
decomposition of IS transitions onto a basis of modes. 
In Section [V] we will discuss the quantitative results of 
the decomposition procedure just described, and finally 
in Section I VII we will summarize our results and frame 
them in a broader context. 



II. MODEL DETAILS 

We study the two dimensional "65:35" Kob- Andersen 
binary Lennard- Jones mixture^ 7 - In brief, it is character- 
ized by the parameters oab = 0.8ctaa, obb = 0.88ctaa, 
(,ab = 1-5caa, and esB — with 65% of the 

particles being of type A. All particles have unit mass. 
Time scales are reported in Lennard- Jones units with 
t = y/ ma AA I €aa and we fix the number density at 
p = 1.2. All interactions were shifted and cut off at 
r = 2.5cTjj. This model has been used previously be- 
cause it resists crystallization in two dimensions better 
than the standard 80:20 variant . 37 ' 38 To our knowledge, 
structural and dynamical quantities for this system have 
not been previously reported. Here, for completeness, we 
briefly discuss these properties for N = 1000. In Fig. Q] 
we present the radial distribution functions for a series 
of temperatures that span the high and supercooled tem- 
perature regimes. In Fig. [2ja) we show the self-part of 
the intermediate scattering function for the A-type par- 
ticles, Fg(k = 6.28, t), as well as the a relaxation times 
as defined by F s (k = 6.28, t = r a ) — 1/e. These val- 




FIG. 2. (a) Self-intermediate scattering function for A type 
particles from NVE molecular dynamics at the temperatures 
listed in Table [I] appearing from left to right. Inset: Alpha- 
relaxation times for both A and B type particles, with data 
given in Table [T] (b) F^(k = 6.28, t) with the horizontal axis 
scaled by the alpha relaxation time r^. The curves collapse 
in the /3-regime as the temperature is lowered. The highest 
three temperatures are not shown. The dotted line shows the 
value 1/e. 
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TABLE I. Alpha relaxation times r Q for the 2D, 65:35 Kob- 
Andersen system, with N = 1000. Reported for A and B type 
particles, as well as for all particles. Also listed for reference 
are integration time-steps used both in annealing the config- 
urations and in generating inherent structure trajectories. 



ues are also reported in Table HI Finally, in Fig. [2jb) we 
also show that time-temperature superposition appears 
to hold in the /3-regime in a manner comparable to the 
3D, 80:20 variant of Ref. M- 



III. INHERENT DYNAMICS 

The motion of particles in a supercooled liquid is ex- 
tremely complex, and the ability of a single particle to 
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FIG. 3. (a) Inherent structure energies (per particle) at T = 0.7 (above) and T — 0.5 (below), plotted as a function of time 
scaled by the alpha relaxation time r a at that temperature. Transitions occur between extended plateaus at both temperatures; 
they occur less frequently on an absolute time scale at the lower temperature but more frequently relative to the alpha relaxation 
time. The dashed box shows the area plotted in (b). (b) Above is a magnified look at the transitions at T — 0.5, with energies per 
particle shown on the left. Below are the maximum squared displacements between subsequent inherent structures, with values 
shown on the right. Note every energy transition has a structural difference as measured by maximum squared displacement, 
and there is an additional transition at 10t/r a ~ 10.8 not visible in the energy trace. 



change its current position is dependent on the interplay 
of fluctuations on many length scales. The notion of in- 
herent structures (IS) has been developed as an aid to 
simplify the description of both structural and dynami- 
cal processes in supercooled liquids42r— In the IS frame- 
work, the motion of the system as a whole is decomposed 
into transitions between local minima on the global po- 
tential energy surface (the "inherent structures" ) and the 
vibrations around these configurations. A temporal se- 
ries of IS generated by quenching a trajectory from a 
molecular dynamics simulation is termed an IS trajec- 
tory. Studies of these trajectories have provided interest- 
ing insights into the true dynamics in glassy systems, for 
example, revealing the presence and nature of string-like 
cooperative motion42r— 

For this study, we will require IS trajectories at a se- 
ries of decreasing temperatures. Specifically, configura- 
tions of 250 particles each were generated using NVT 
dynamics in LAMMPS by annealing through the series 
of temperatures in Table HI simulating for over 100r Q 
(as calculated for N = 250) at the lower temperature for 
each subsequent annealing step 4^ Each resulting config- 
uration was then simulated for an additional 100r Q , and 
F s (k,i) was calculated for that entire trajectory, as well 
as for the first and second halves. At the lower tempera- 
tures, signs of aging were obvious in the majority of these 
simulations; only those configurations resulting in a sim- 
ulation that showed little or no difference in dynamics 
between the first and second half of the trajectory were 
used for subsequent analysis. At all temperatures, we ob- 
tained at least 20 independent configurations for future 
study. The integration time steps used throughout this 
work are reported in Table |TJ 

For each of these independent configurations, IS trajec- 
tories were generated by quenching configurations from 
molecular dynamics simulations using steepest descent 



minimization. Our system is quite small, rendering the 
IS meaningful in concept and easy to identify. Specifi- 
cally, trajectories of length 100,000 integration steps were 
generated at each temperature, saving every configura- 
tion, except at T = 0.45 where we chose to run 200,000 
and save every 2 steps as a means of obtaining better 
statistics. We generated IS trajectories using a bisection 
technique similar to that of Ref. 44l A small number of 
configurations evenly spread throughout the trajectory 
were quenched to a local minimum. If two adjacent con- 
figurations differed (see next paragraph for more details) , 
then we bisected that time interval and quenched the 
configuration halfway between the two endpoints. This 
procedure was iterated until a sequence of configurations 
was generated in which all transitions between minima 
occurred at the time scale of a single integration time 
step. 

For every configuration in the new IS trajectory, 
we monitored both the energy and the displace- 
ment from the previous configuration. We calcu- 
lated both the mean square displacement (ARjg(tj)) = 

-ffY^rj=i\ r j(ti) ~ r j(ti-i)\ 2 an d maximum squared 
displacement max{(<5rf) 2 } = maxi<j<jv{|ij(ti) — 
r j{ti-i)\ 2 } where £, is the time index of the i'th inherent 
structure, r^(tj) is the position of particle j in inherent 
structure i, and N is the total number of particles. We 
find that these two quantities are correlated and signifi- 
cant changes in one always occur in the other simultane- 
ously. 

A plot of the IS energy, as displayed in Fig. (Ha) , re- 
veals a series of plateaus as in previous worksi 44 ' 47 " — 
Differences in energy within these plateaus are on the 
order of 10 -6 per particle or smaller. However, as dis- 
cussed by Schr0der, et al, there is the small probability 
that a change in structure occurs between two structures 
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FIG. 4. At each temperature studied, the proportion of in- 
herent structure transitions histogrammed into bins of size 
P(AR I ) — 0.1. At high temperatures there are signifi- 
cantly more "extended" transitions than at low temperatures. 
Dashed lines are a guide to the eye. Note that calculations 
on transitions with P(AR I ) > 0.55 such as those in Figs. [9] 
and [10] may be skewed due to poor statistics, but this does 
not affect any of our conclusions. 



with a vanishingly small difference in energy. 44 Hence, 
as in that previous work, we will identify changes in in- 
herent structure by the difference in positions and not 
by energy. We used max{(Jr 7 ) 2 } for this purpose but 
emphasize that choosing the mean square displacement 
would produce identical results. For the model studied 
here, we find that maximum squared displacements in 
the range of 10 -5 to 10 -2 are very rare. We therefore 
empirically choose a threshold value within this range 
and deem that whenever max{(<5r 7 ) 2 } > 0.001, a tran- 
sition has occurred. The correlation between changes in 
structure and changes in energy is illustrated in Fig.(3Jb). 

It will be most important in future discussion to quan- 
tify the degree of localization of motion in an IS transition 
(or a normal mode, see Sec. IIV|) . The metric we shall use 
is the standard participation ratio. For a vector quantity 
v of length dN such as the displacement field, the vec- 
tor can be organized as a matrix V of size N x d where 
each row 1 < i < N contains a quantity pertaining to 
particle i. With these quantities, the participation ratio 
is defined asJ^ 



P(v) 



JV 



A^(V 4 -V 4 ) 2 



(1) 



We use this definition to calculate the participation 
ratio of transitions P(AR I (U)) at each temperature. 
Nearly all participation ratios fall in a range from to 
0.6, with just a few occurring with values larger than 0.6. 
As a matter of nomenclature only, we will refer to tran- 
sitions with F(AR / (t l )) > 0.3 as "extended", and all 
others as "localized". A histogram is shown in Fig. 2) As 
the temperature is lowered, transitions with low partici- 
pation ratio become much more prevalent i.e. motion is 
increasingly localized. Interestingly, the curves intersect 
for participation ratio w 0.17, although the meaning of 



this coincidence is unclear. An illustration of the types 
of IS transitions is shown in Fig. [5] 



IV. NORMAL MODES 

Suppose a transition in a given IS trajectory is found 
to occur at time i,. We wish to study the relationship 
between the normal modes of configuration TLjs(U) and 
the transition Rjs(ij) — > R/s(ij + i). To do this, we con- 
struct the Hessian or dynamical matrix for configuration 
R/s(^) is given by the matrix of second derivatives of 
the potential 



H(ii 



d 2 V(H IS (U)) 
drjdrl 



(2) 



Normal modes are the eigenvectors obtained by diago- 
nalization of this matrix. It is well known that the normal 
modes form a complete orthonormal basis for the motion 
of this configuration. For configuration i the eigenvectors 
can be ordered by eigenvalue and will be denoted {e*}. 
For a d-dimensional configuration of N particles, e* will 
be a vector of length dN and there will be dN eigenvec- 
tors, of which d will correspond to translational modes 
with eigenvalue 0. The Hessian is real and symmetric and 
hence the eigenvalues and eigenvectors are real. The den- 
sity of states (DOS) D(u>) is a histogram of frequencies 
obtained from the eigenvalue spectrum, and is plotted for 
our system at a series of temperatures in Fig. [5fa) . The 
same histogram is shown with the DOS divided by a fac- 
tor of oj in Fig. OJb) so that one can identify the so-called 
"boson peak" , where the number of low frequency modes 
exceeds the Debye prediction of 

The displacement vector for a transition is given by 
|AR 7 (^)) = \His(ti)) — |R/s(f,_i)), and characterizes 
the inherent structure transition. For convenience, we 
will consider instead the normalized displacement vector 
defined as |AR 7 (i 2 )) = AR 7 ^)/^ (AR 7 (ti)|AR 7 (^)). 
This can be projected onto a normal mode basis as in 
Eq.l 



dN 

|AR 7 (t 4 ))=E c "Wl e ") 



(3) 



Note that we label the eigenvectors with configuration 
index k to emphasize that this need not be the eigenval- 
ues corresponding to time fj. The coefficients dj(i) above 
are obtained through the orthonormality condition for 
the basis 



dN 

(e^|AR 7 (t 4 ))=^cf( l )(e J fe |ef)= S fe ( l ). 
i=i 



(4) 



Since |AR 7 (ti)) is normalized, its squared norm can 
be written as 
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3 = 1 



(5) 
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(a) 



FIG. 5. (a) Density of states for the 2D Kob- Andersen model 
with N = 250 at the series of temperatures shown. Curves are 
shifted up by 0.002 from the temperature below for clarity, the 
lowest temperature being at the bottom, (b) The same as (a) 
but the density of states has been divided by frequency, and 
scaled up by a factor of 10 for clarity, (c) Participation ratio 
as a function of frequency for the same series of temperatures 
as in (a). Curves are shifted up by 0.1 from the temperature 
below for clarity. Quantities in (a) and (b) are binned to a 
resolution of 0.2r _1 . 



We will refer to |c^'(i)| 2 as the contribution of mode 
j of basis k to transition i. These contributions can be 
ordered from highest to lowest in numerical value. 

In order to get an idea of how much each mode con- 
tributes to a transition, we consider the quantity Ni/ 2 , 
which is defined to be the smallest number of modes 
necessary to reproduce 50% of a given transition. We 
will also report F 1 / 2 , the fraction of all modes necessary 
to reproduce 50% of the motion, which has been used 
previously . 24 ' 36 In general we can also define the quan- 
tities Nf and Ff where / is a fraction other than 0.5. 
To provide a unique definition of these quantities, we 
take {of(i)} as a monotonically decreasing sequence of 
ordered contributions to transition i from basis k and 
then generate the cumulative series sj^(i) = SIHi °/ fc (*)- 
Then — 



Nf[i) = min{m|«*,(i) > /} 



(G) 



and Ff(i) = N^{i)/{Nd). 

The definition of participation ratio given in Eq.[T]also 
applies to the normal modes obtained via diagonalization 
of the Hessian. By convention, the participation ratio 
averaged over modes with a given frequency is referred 
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FIG. 6. Inherent structure transitions for the 2D Kob- 
Andersen model with N = 250 at T — 0.5. Axis labels are 
simulation coordinates, (a) A "localized" transition with par- 
ticipation ratio P(AR/) = 0.04 and ((AR 7 ) 2 ) = 0.06. Par- 
ticle displacements are shown with black arrows and a pro- 
jection of the top 37 highest contributing modes, reproducing 
50% of the transition is shown with blue arrows. Vectors 
are multiplied by a factor of two for clarity, (b) The same 
transition with the natural logarithm of square particle dis- 
placements ln((<5r 7 ) 2 ) represented from red at —12.5 to blue 
at —0.7. (c) Same as (a) but for an "extended" transition with 
P(AR J ) = 0.5 and {(AR 7 ) 2 ) = 0.2. Here only one mode is 
needed to project onto 60% of the transition, (d) The same 
transition as (c), with ln((<5r 7 ) 2 ) ranging from —9.43 in red 
to -1.48 in blue. 



to as P u . This quantity, which will play a role in our 
future discussion, is shown in Fig. [3Jc) for a wide range 
of temperatures. 

In Fig. [S] we illustrate the concepts presented in this 
and the preceding section. Fig. [Ha) and (b) show an 
IS transition with low participation ratio, in which most 
of the contribution to the motion is localized to a few 
particles. In Fig. (SJc) and (d) we show an IS transi- 
tion with high participation ratio, in which the motion is 
much more evenly spread throughout the system. In (a) 
and (c) we show the result of projecting the motion of 
an IS transition onto normal modes of the initial IS. The 
black arrows show the motion of the particles in the sys- 
tem during an IS transition, and the blue arrows show 
the sum of the normal modes that contribute to N 1 / 2 
weighted by their coefficients as in Eq. [3] We see that 
for both localized and extended transitions, the degree 
of similarity is striking. In the single example of an "ex- 
tended" transition, a single normal mode projects onto 
60% of the actual displacement vector. We illustrate the 
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opposite case with a particularly localized event, where 
just a few particles rearrange, and where a relatively large 
number of modes are needed to reproduce half of the mo- 
tion compared to the average as discussed in Sec|V] In 
what follows, we use the metrics N 1 / 2 and F 1 / 2 as a quan- 
titative measure of how faithfully a given transition may 
be represented by a small set of normal modes. 



V. RESULTS 

We begin this section with a discussion of our main 
result. We examine the temperature dependence of the 
quantity (Ni/ 2 ) as defined in Section HVl (where brack- 
ets denote average over many transitions) following Brito 
and Wyart, who showed that this quantity decreased sys- 
tematically in their system of hard spheres as the jam- 
ming transition was approached. 36 Analogously, we plot 
in Fig. [JJ the identical quantity for our thermal system 
as a function of inverse temperature. It is clear that as 
our system becomes increasingly supercooled the num- 
ber of normal modes needed to approximate the corre- 
lated atomic motion occurring in an inherent structure 
transition decreases. The degree to which this is true is 
independent of whether the pre- or post-transition inher- 
ent structure is used to produce the normal mode basis. 
This fact, which is required if we are to interpret the 
mode basis as a set of directions along which reversible 
motion occurs on a potential energy landscape, is not en- 
tirely trivial as the modes of the two configurations are 
distinct 

To gain some insight into the magnitude of (N 1 / 2 ), as 
well as the significance of the trend seen in Fig. [JJ we 
also reconstruct the dynamical activity exhibited during 
inherent structure transitions with a random set of modes 
obtained from distant inherent structures sampled at the 
same temperature as the transitions under consideration. 
The inset of Fig. [Jj shows that the value of (iVm) ob- 
tained from this procedure exceeds that obtained with 
modes associated with the actual transition by approxi- 
mately one order of magnitude at the temperatures con- 
sidered in our simulations. Further, as temperature is 
lowered, the average number of random modes needed to 
approximately reproduce the motion observed in inherent 
structure transitions remains constant, in stark contrast 
to the behavior seen when (Ni/ 2 ) is calculated with the 
modes associated with the IS transitions. In Fig. [5] we 
show that behavior illustrated in Fig. [JJ does not depend 
on the choice of the threshold value / used to define the 
degree of overlap between the modes and the true par- 
ticle motion. The behavior illustrated in Fig. [Jj suggests 
that motion in mildly supercooled liquids proceeds along 
a successively smaller number of "soft" mode directions, 
a phenomena noted previously using a different definition 
of normal modes as hard spheres approach the jamming 
transition^ 6 - 

The behavior illustrated in Figs. [JK] is, upon first 
consideration, contradictory with the result presented 
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FIG. 7. Average number of modes (left axis) or fraction of 
total modes (right axis) necessary to reproduce at least 50% 
of total particle motion in inherent structure transitions at 
each temperature studied. Blue squares are for transitions in 
the "forward" direction and red circles in the "backwards" 
direction as described in the main text. Error bars shown 
are jackknife standard errors.— Inset: Symbols are the val- 
ues of {Ni/ 2 ) and {-F1/2) where 2000 inherent structures at 
each temperature are projected onto the modes from each 
of twenty-five random configurations from an independent IS 
trajectory. Note the difference in magnitude between the ver- 
tical axis of the inset and the main figure. 
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FIG. 8. Average number and fraction of total modes neces- 
sary to reproduce a fraction / of inherent structure transitions 
from / = 0.3 to / = 0.8. A dotted line shows / = 1/2, the 
value focused on for most of this work. A slice through this 
line gives the data shown in Fig. [7] (Nf) decreases with tem- 
perature for all values of /. Jackknife error bars are approxi- 
mately the size of the symbols.— Dashed lines are a guide to 
the eye. 



in Fig. [Jj In particular, Fig. [4] shows that the motion 
that occurs during IS transitions becomes increasingly 
localized as temperature is decreased. In Fig. |B]we illus- 
trate that, at least in one particular transition, it takes 
more modes to construct (N1/2) for the case of a local- 
ized transition than for a delocalized one. This fact can 
be quantified, and we show in Fig. |H]that this is true on 
average at any temperature. Thus it is surprising that 
(-W1/2) decreases as temperature is lowered. The resolu- 
tion of this conflict lies in the fact that, while for a given 
temperature(iV 1 /2) is larger in localized cases, the curves 
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FIG. 9. {N1/2) and {F1/2) for inherent structures as a 
function of participation ratio of IS transitions. Bins of size 
P(AR J ) = 0.1 are used. "Extended" transitions (measured 
by P(AR 7 ) require fewer modes to reproduce their transi- 
tions. (iVi/2) and (-F1/2) decrease uniformly for all values of 
P(AR 7 ) as temperature is lowered. Jackknife error bars are 
again approximately the size of the symbols and dashed lines 
a guide to the eye.— Data for P(AR 7 ) > 0.55 are not shown 
due to poor statistics, as seen in Fig. [4] 
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of (N1/2) uniformly decrease for all types of transitions 
as temperature is lowered. Clearly either something in 
the structure of the modes or the type of particle motion 
occurring in IS transitions (or both) changes as tempera- 
ture is lowered such that the overlap between motion and 
modes increases in a uniform fashion independent of the 
degree of localization of the transition. Any changes in 
the mode structure are evidently subtle, at least in terms 
of the common metrics shown in Fig. [5] In particular, 
since the temperature dependence of both the density 
of states and density of participation ratios of IS normal 
modes is very weak, differences must be found at the level 
of the small number of modes that actually contribute to 
(-W1/2) & t each temperature. 



FIG. 10. (a) Average frequency of modes contributing to 
{N\/2) as a function of IS transition participation ratio his- 
togrammed into bins of size P(AR 7 ) = 0.1. "Extended" 
transitions are made up of lower frequency modes. The fre- 
quency of contributing modes decreases with temperature for 
all sized transitions. Jackknife errors are smaller than the 
points shownjSi Dashed lines are a guide to the eye. (b) 
Frequency distribution of modes used to compute (u>) for 
three temperatures T = {0.9,0.6,0.45}, and two values of 
P(AR J ) = {0.05,0.45} (referred to as P here for conve- 
nience). Other temperatures and P values interpolate be- 
tween these curves. Temperatures are shifted to the right from 
the temperature above by ui — 10 and curves for P = 0.45 are 
dotted. We see that curves are more sharply peaked for higher 
P and the peak height increases as temperature is lowered. 



While the overall features of the various densities of 
states are largely insensitive to the lowering of the tem- 
perature, we can show that the properties of the modes 
that actually contribute to(iV 1 / 2 ) do have properties that 
systematically change as temperature is varied. The con- 
tributing modes are essentially confined to the low fre- 
quency band in the neighborhood of the "boson" peak 
(with frequencies of 10r _1 or smaller) at all temperatures 
studied. Within this band, there is a systematic shift of 
the frequencies of contributing modes to lower values at 
lower temperatures, as illustrated in Fig. 1101 Further- 
more, for localized transitions the fraction of localized 
harmonic modes that contribute to {N 1 / 2 ) systematically 
increases as temperature is lowered, as shown in Fig. 1111 
Thus, as temperature is decreased the modes that con- 
tribute to (N1/2) increasingly take on the character of 
the actual motion of particles during the IS transition. 



VI. CONCLUSION 

The results of the previous section clearly show that 
dynamical heterogeneity in supercooled liquids, as ex- 
hibited by the motion of particles during IS trajectories, 
may be cleanly decomposed into a small number of har- 
monic normal modes associated with the inherent struc- 
tures themselves. The fraction of Hessian eigenvectors 
needed to reproduce an average IS transition as quanti- 
fied by (N1/2) is on the order of two percent of the total 
number of modes and decreases as the liquid is increas- 
ingly supercooled. This decrease in the number of soft 
mode directions that the system follows takes place even 
though IS transitions become more localized. Further, 
we have demonstrated that this property is a nontrivial 
feature of the actual pre- and post-IS mode structure and 
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FIG. 11. (a) Average number of modes which have a partic- 
ipation ratio < 0.4 and which contribute more than 5% 
to transitions of a given type at three different temperatures, 
(b) For the same temperatures, average total contribution to 
a single transition of a given type from modes with < 0.4. 
Error bars are jackknife standard errors.— Modes with low 
participation ratios contribute more to "localized" transitions 
(P(AR 7 ) < 0.3) at low temperatures than at high tempera- 
tures. 



cannot be reproduced by modes obtained from indepen- 
dent inherent structures. 

This study was clearly inspired by the work of Brito 
and Wyart, who found similar behavior in hard-spheres 
as the packing fraction is varied and the jamming tran- 
sition is approached ! 24 ' 36 Our work shows that this re- 
lationship between modes and dynamical heterogeneity 
is similar in thermal systems where inverse temperature 
plays the same role as density. Interestingly, we find that 
for an analogous range of relaxation times, even fewer 
modes are needed in our study to reproduce the quan- 
tity (Ni/ 2 ) than needed in the hard-sphere case. This 
is somewhat surprising because the statically defined IS 
normal modes are, in a sense, cruder than the dynam- 
ically defined modes of Brito and Wyart which are as- 
sociated with the free energy (as opposed to potential 
energy) of the system. It is still unclear if the results 
presented here can be taken as a feature of similarity be- 
tween the jamming and glass transitions. In particular, 
given the importance of marginal stability in jamming, 
the depletion of mode directions as a basis for motion 
in jamming is natural. On the other hand the role of 
such mechanical consideration in thermal glass-forming 
systems is totally unclear. One interesting possibility 



hinted at by the lowest temperature data point in Fig. [7] 
is that the temperature dependence of (N 1 / 2 ) plateaus 
at low temperature, perhaps close to the mode-coupling 
temperature of the system i 21 ' 52 ' 53 While the quality of 
statistics we have is not good enough to draw such a 
conclusion, this possibility is worthy of future scrutiny. 
In particular, the jamming transition is accompanied by 
a divergence of the length scale associated with mode 
polarization and hence the directions associated with re- 
laxation vanish at the transition* 1 ^ In thermal systems, 
this behavior is absent or averted due to other processes, 
leading to a possible saturation of the number of soft 
directions that overlap with the true relaxation of the 
system. 

Another avenue worthy of future study concerns the 
degree to which the results presented here are an indi- 
cation of the primacy of soft modes in the relaxation 
of supercooled liquids. In an interesting study, Ashton 
and Garrahan constructed models with mode structure 
by taking a standard lattice gas and connecting particles 
via springs with random force constants.— They showed 
that, within models where soft- modes are grafted onto 
kinetically constrained lattice gasses, the soft mode po- 
larization will indeed congregate near vacancies but do 
not govern the motion of the defects, which instead move 
via local facilitated rules. Thus, in this class of models, 
soft modes are passive markers of defect structure that 
do not themselves participate in relaxation. It would 
be most interesting to repeat the exercise of decompos- 
ing the inherent structures transitions in kinetically con- 
strained models onto the normal modes of the Ashton- 
Garrahan models to test if the picture of modes as pas- 
sive markers is compatible with the behavior exhibited 
in Fig. [7] and the previous work of Brito and Wyart. We 
reserve this investigation for a future study. 
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